Coherent control of localization, entanglement, and state superpositions in a double 

quantum dot with two electrons 
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We have recently proposed a quantum control method based on the knowledge of the energy 
spectrum as a function of an external control parameter [Phys. Rev. Lett. 99, 036806 (2007)]. So 
far, our method has been applied to connect the ground state to target states that were in all 
cases energy eigenstates. In this paper we extend that method in order to obtain more general 
target states, working, for concreteness, with a system of two interacting electrons confined in 
semiconductor double quantum wells. Namely, we have shown that the same basic method can 
be employed to obtain localization, entanglement, and general superpositions of eigenstates of the 
system. 

PACS numbers: 7,3.63.-b, 78.67.Hc 
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I. INTRODUCTION 

Quantum control is an area of research of great current 
interest with a vast potential for applications in quantum 
information technologies lii^ The basic problem of quan- 
tum control consists of driving externally a quantum sys- 
tem in order to take it from an initial state to a target 
final state. Quantum control ideas and methods are tra- 
ditionally widely applied in magnetic resonanceS^i and 
quantum chemistry,^ and they are advancing rapidly in 
the area of soHd state nanostructures.i^>^>^ 

In a series of recent publication a^°i^^d^ we have pre- 
sented a method of quantum control based on the knowl- 
edge of the energy spectrum as a function of a suit- 
able single control parameter. This method is useful 
provided that the transitions between neighboring lev- 
els are well described by the Landau-Zener model.— The 
latter condition has allowed us to successfully navigate 
the spectrum with a combination of diabatic and adia- 
batic changes of the control parameter. Other authors 
have also employed the navigation of the energy spec- 
trum as a coherent quantum control tool, especially in 
the field of atomic and molecular systems controlled by 
lasers.— li^ii^iiiii^ We have applied our control method 
to two different physical systems with remarkable suc- 
cess. The first one, which we will study in this paper, is 
a nanostructured semiconductor system consisting of two 
interacting electrons confined in quasi-one-dimensional 
quantum dots,— »ii>i^ and the second one was the LiCN 
moleculepi^ In both systems we have been able to con- 
nect distant eigenstates through long and complex paths 
in the energy spectrum with very high probability. 

In this paper we take the control method farther: by 
allowing not only diabatic and adiabatic transitions but 
also intermediate velocities, we can arrive at more com- 
plex final (target) states. With this generalized naviga- 
tion method we are able to achieve the creation of lin- 
ear superpositions of adiabatic states starting from the 
ground state. This opens a rich menu of possibilities, like 
the creation of Bell states and other types of entangled 
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Figure 1: Schematic plot of the building block of our con- 
trol method. An avoided crossing and a schematic drawing of 
the diabatic (fast), adiabatic (slow), and intermediate ways 
to cross it. (a) Long (short) arrows represent fast (slow) vari- 
ations of the control parameter; (b) Medium arrows represent 
the intermediate velocity. 



states. 

The paper is organized as follows: in Section [Ul we de- 
scribe the main ingredients of our control method. In 
Section IIIII we present again the two-electron nanostruc- 
ture. Section [IV] is devoted to the results of the gener- 
alized control method. Among the results presented, we 
show how to apply our method to generate Bell states 
and coherent superpositions of energy eigenstates. 



II. THE CONTROL METHOD: NAVIGATING 
THE HILBERT SPACE 

Let us first review the basic ideas of our method in the 
simplest possible system, i.e. a parameter-dependent two 
level system. Let us assume that these two levels have an 
avoided crossing as shown in Fig.[TJ At the avoided cross- 
ing, the two energy levels approach each other and the as- 
sociated eigenstates exchange their characteristics as the 
external parameter sweeps through the crossing. If the 
avoided crossing is traversed very slowly (adiabatic path 
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Figure 2: Confining double well potential in the longitudinal 
direction of the coupled quantum dot structure. The external 
electric field is -B = (solid lines) and E = 12 kV/cm (dashed 
lines). 



in Fig.[T][a)), the adiabatic theorem guarantees that one 
will stay in the initial adiabatic level, but with the para- 
doxical consequence that the final state will have been 
exchanged with the other diabatic state. On the other 
hand, if the avoided crossing is traversed very quickly 
(diabatic path in Fig. [1]), the characteristics of the state 
are preserved. Clearly, these two limiting possibilities 
suggest a very simple control method. The quantitative 
meaning of slow and fast transitions is given by the the- 
ory of Landau-Zener transitionsj^i-^^ 

This basic idea has allowed us to perform compHcated 
control tasks. Namely, we can travel through the energy 
spectrum using the avoided crossings as sources of effi- 
cient controllable choices between two adiabatic states. 
In this way we can connect distant adiabatic states pro- 
vided that there exists a path going from one to the other 
via jumps at avoided crossing and adiabatic evolutions in 
the absence of crossings. Examples of the applic ation of 
this method have been presented in Refs. |lC)llll| . 

In this paper we go one step further by using not 
only diabatic and adiabatic transitions but also transi- 
tions with intermediate speeds. This type of transition 
gives final states (on exit of the avoided crossing) which 
are linear combinations of the two adiabatic states (Fig. 
[ll[b)). The combination of several crossings using inter- 
mediate velocities allows us to access a great variety of 
final states. Thus, in this work we generalize the control 
method proposed earlier and present several applications 
which show the power of the improved control technique. 



III. THE DOUBLE QUANTUM DOT SYSTEM 

In order to explore on a concrete system the idea of us- 
ing intermediate speeds to traverse avoided crossings we 
will study here a system of two interacting electrons in- 
side a quasi-one-dimensional double quantum dot which 
was used in our previous works. The system is subject 
to an external, uniform electric field, whose amplitude is 
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Figure 3: Main panel: The energy spectrum of the two inter- 
acting electrons confined in a quasi-one-dimensional double- 
well semiconductor nanostructure as a function of an external 
uniform electric field. The first 31 energy levels are shown. 
Side panels: Spatial wave functions 4>i{E,zi,Z2) correspond- 
ing to labels (a) to (h) of the main panel. States (a) and (b) 
have one electron in each well, wave functions (c) and (d) are 
localized in the left well, and wave functions (e) and (f) have 
both electrons in the right well. States (a) to (f) are far from 
avoided crossings and therefore have well-defined localization 
properties. This is not the case for eigenstates (g) and (h), 
which are at avoided crossings. 



used as the control parameter. The confinement of the 
two electrons is very narrow on two dimensions, which we 
denote as x and y, and the double well profile is defined 
along the remaining, longitudinal coordinate, z. The ef- 
fective Hamiltonian of the two electrons reads 



H 



1^ 1^2+1^2 +V{ZI)+V{Z,) 



2m \dzf 
+Vci\zi ~ Z2\) - e{zi + Z2)E{t), 



dzl 



(1) 



where m is the effective electron mass in the semiconduc- 
tor material, Vc is the effective Coulomb interaction be- 
tween the electronSj^ft is the confining potential (see 
Fig. [21), and E{t) is the external time-dependent electric 
field. We choose as confining potential a double quantum 
well with well width of 28 nm, interwell barrier of 4 nm, 
and 220 meV deep (a typical depth for a GaAs-AlGaAs 
quantum well) . We normally assume that the initial state 
is the ground state, which is a singlet. Since the Hamil- 
tonian is spin independent, the total spin is conserved 
and the spatial wave function remains symmetric under 
particle exchange at all times. 

Let us review some characteristics of the energy spec- 
trum of the two electrons and of the eigenstates which 
will guide our control strategies. First we consider the 
spectrum as a function of the control parameter (the ex- 
ternal electric field). This spectrum is plotted in Fig. 
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[3] with selected eigenstates. The energies ei{E) and the 
eigenstates (j)i{E, zi,Z2) are obtained by numerical diago- 
nalization. For this calculation, we have used as basis set 
the symmetric combinations of the twelve lowest single- 
particle eigenfunctions. Thus, our basis set of the two- 
particle Hilbert space has 12*(12+l)/2=78 states. An 
important characteristic of this system in the lower part 
of its spectrum is that the states have fairly well-defined 
localization properties. The three possible types of local- 
ization (both electrons in the left dot, in the right dot, or 
electrons in different dots) are associated with the three 
types of slopes seen in the spectrum (positive, negative, 
or almost zero, respectively). In Fig. [3] we show some 
states to illustrate the mentioned locaHzation character- 
istics. In states a and b the electrons are in different dots, 
in states c and d both electrons are in the left dot, and 
in states e and f both electrons are in the right dot. Of 
course, at the avoided crossing the states get mixed and 
these well-defined properties are lost (states (g) and (h)). 



IV. RESULTS 
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Figure 4: Probability Pll as a function of the time-dependent 
electric field i?(t). The velocities of the electric field are: (i) 4, 
(ii) 1, (iii) 0.4, (iv) 0.2, and (v) 0.04 (kV/cm)/ps. With circles 
we represent Pll of the adiabatic ground state as a function 
of the electric field (see text for details.) Inset: closer view 
of the first avoided crossing involved in these processes. The 
arrows indicate the adiabatic path. 



This section is devoted to illustrate the power of the 
generalized control method in which the velocities to 
cross the avoided crossings are not restricted to produce 
diabatic and adiabatic transitions. In other words, we 
will show the possibilities opened by the ability to allow 
intermediate velocities. Due to the potential for appli- 
cations of our method, we consider it best to illustrate 
its power through two important examples of coherent 
control. 



A. Localization 

An application of quantum control that has been ex- 
tensively explored in the recent literature is the locaHza- 
tion of one or two electrons in a double well potential. 
The general idea is to start the coherent evolution in the 
ground state, which in the case of two electrons is de- 
localized due to the PauH and Coulomb repulsions, and 
end up in a state in which both electrons are in the same 
well. In order to describe the degree of localization of the 
two electrons we define the probability for both electrons 
to be in the left well, Pll- 

/ dzidz2|V(^l,22,t)|2, (2) 

-OQ J ~OQ 

where ip{zi,Z2,t) is the evolving wave functions of the 
two electrons. 

An easy way to localize both electrons in the left dot 
using our control method starting from the ground state, 
consists of varying the electric field adiabatically in or- 
der to pass the first avoided crossing between the first 
two levels located at E = 4.81kV/cm. This process is 
shown in the inset of Fig. [H (The numerical solution of 



the time-dependent Schrodinger equation was obtained 
using the usual fourth-order Runge-Kutta method with 
a time step of 0.05 fs.) The electric field varies linearly 
with time and we display the evolution of Pll for several 
velocities. We start at the ground state, and accordingly 
Pll — at t = 0. For a high velocity, we expect to 
cross diabatically the avoided crossing and the locaHza- 
tion properties will not change considerably. This case is 
seen in curve (i) {E = 4 (kV/cm)/ps.) in Fig. |4l As the 
velocity is decreased (curves (ii) to (v)) we approach a 
final state that is highly localized in the left dot. There 
is a maximum value of the probability Pll that can be 
obtained with this method, given by the probability Pll 
of the adiabatic ground state after the avoided crossing. 
The probability of the adiabatic ground state is plotted 
with open circles in Fig. |4l After the crossing it is a 
slowly rising function in the plotted range and reaches 
the value of 0.998 at electric field E = 9kV/cm. We see 
that for the velocity of curve (v) {E — 0.04 (kV/cm)/ps) 
the evolving probability Pll(0 follows tighlty the open 
circles and near = 9 kV/ cm presents a small oscillation 
bound between 0.995 and 0.9975. That is, the maximum 
value of 0.998 is approached within a 0.3 percent. 

Another simple way to obtain localization in this sys- 
tem is what we may caH the sudden-switch method, 
which uses step-wise constant fields;^ We wish to com- 
pare the effectiveness of both methods to obtain localiza- 
tion. First, let us summarize the basic procedure of the 
sudden-switch method (See Fig. [5]). Starting from the 
ground state one applies two successive steps of constant 
electric field. In the first step one goes from zero field 
to the field corresponding to the first avoided crossing, 
i.e. E = 4.81kV/cm in our case. While the new field is 
on, the probability Pll oscillates with the frequency cor- 
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Figure 5: Localization with the sudden-switch method. 
While the electric field stays at the avoided-crossing value 
{E = 4.81 kV/cm) the probability PlL oscillates with the fre- 
quency corresponding to the energy gap (dashed lines). The 
method consists of increasing suddenly the electric field when 
PlL is maximal, freezing the localization on its highest value 
(full lines). Upper panel: Probability that both electrons are 
in the left dot. Lower panel: Step-wise constant electric field." 



responding to the energy splitting at the avoided cross- 
ing. This oscillation occurs because the initial state is 
no longer eigenstate of the Hamiltonian at the avoided 
crossing. In fact, it is a 50%-50% linear superposition 
of the two adiabatic states (eigenstates) at the avoided 
crossing. When the probability Pll reaches a maximum, 
that means that the current state is the other diabatic 
state. At this time one switches again the field to leave 
the avoided crossing. Far from the avoided crossing the 
diabatic states are very close to the eigenstates of the 
Hamiltonian, and therefore, are approximately station- 
ary. Thus, the Pll remains at the highest value attained 
during the oscillations at the avoided crossing. With this 
method, a localization Pll of up to 93% can be achieved, 
as can be seen in Fig.[5j While the time scales involved in 
both methods are the same, in comparison, our method 
has three advantages: (i) a higher degree of localization 
can be obtained; (ii) the sudden-switch method requires 
a fine adjustment of the timing not needed in our case; 
(iii) our method is more powerful in the sense that can 
be used to navigate in the spectrum and connect distant 
statesjAft 

Let us now consider states with a different and in a 
sense more complex kind of localization. As mentioned 
earlier, we can find in the spectrum three types of local- 
ized states, which we can denote in the following way: 
\RR), \LL), and \RL), which have, respectively, both 
electrons in the right and left dots, and one electron 
in each dot. (Of course, these states cannot be consid- 
ered to be product states of single-particle orbitals, be- 



cause quantum correlations are generally present in all of 
them.—) We usually refer to the first two types as local- 
ized states, and to the third one as delocahzed. However, 
a linear superposition of the first two types can also be 
considered as localized, in the sense that one knows that 
both electrons would be found in the same dot if a mea- 
surement were performed. In Fig.[6]we show some control 
paths that may be followed in the spectrum to reach a 
superposition of \RR) and \LL). Linear superpositions 
of the form {\RR) ± \ LL))/ V2 are always available at the 
center of avoided crossings of states with RR and LL lo- 
caHzation. In Fig. [6] (a) and (b) we show how to go from 
the ground state to the Bell- type states ( | RR) -\-\LL))/ y/2 
and {\RR) — |LP))/V2, respectively, using only diabatic 
and adiabatic transitions. We remark that these states 
are energy eigenstates, and thus do not evolve further if 
the electric field is fixed. If the restriction of using only 
diabatic and adiabatic transitions is relaxed we can tra- 
verse the avoided crossing of the two types of states with 
an intermediate speed thus obtaining a combination of 
the diabatic states after the crossing 
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\RR) +be 



-iELLt/h 



\LL). 



(3) 



The values of |a| and |6| can be tuned by choosing the 
appropriate velocity. We note that the effects of using 
intermediate velocities was illustrated previously (Fig.|4]) 
in the context of searching for a LL localized state. 

We now show in Fig. [7] how the path (a) of Fig. [6] is 
obtained by varying the electric field appropriately. In 
Fig. [7] (a) we show the electric field and in (b) we show 
the probabilities as a function of time (semi- log plot). To 
arrive at the desired state, we must cross two avoided- 
crossing adiabatically, one sX E — — 4.81kV/cm and en- 
ergy — — 437meV, and the other one aX E — 3.02kV/cm 
and energy — —402 meV. The first avoided crossing in- 
volves a RL and a RR state, and the second one, RR 
and LL states. The energy gap of these avoided cross- 
ings are very different, and therefore the time required 
to cross them adiabatically is also very different. This 
is way the log scale for the time axis is used in Fig. [7l 
We see that once the electric field is left constant the 
probabilities Pr_r and Pll oscillate between 0.465 and 
0.495 while the overlap of the final state with the target 
one es equal to 0.98, showing that we have arrived at the 
desired state. The other two paths shown in Fig. [6] can 
be realized in a similar way, and with the same degree of 
success. 



B. Coherent superpositions 

In the previous section we used our method to con- 
trol the localization properties of the target state. Now 
we illustrate the fiexibility of our method by obtaining 
coherent superpositions of several adiabatic eigenstates. 
For example, we may wish to obtain a linear superposi- 
tion of the three types of localization having each of them 
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Figure 6: Control paths for two electrons leading to target 
states which are superpositions of states with RR and LL lo- 
calization. The lengths of the arrows indicates the velocity of 
the transitions. The target states are: (a) {\RR) + \LL))/^/2- 
(b) {\RR) - \LL))/V2; (c) state given in Eq. ©. 



the same weight. That is, we seek a state of the form 



1 



\^) = {a\RR) 



b\LL) + c\RL)) , 



(4) 



with \a\ = \b\ = |c| = 1. 

In Fig. [8] the control path to such a state is shown. 
In Fig. [9] we show the overlap of the evolving state with 
the first six adiabatic states and the electric field (the 
control parameter) as a function of time. To appreci- 
ate the evolution in greater detail we plot the evolving 
state at chosen times in Fig. [101 The initial state is the 
ground state (see Fig.[TOl state at time ti =0), and the 
final state is a superposition of the states 1,5, and 6 at 
the electric field E = 3.04kV/cm (see Fig. [lOl state at 
time ie — 1000 ps). We start the evolution by increas- 
ing slowly the electric field (blue arrows in Fig. [8]), so 
that the state remains at the ground state (see Fig. [TOl 
state at time t2), then we accelerate to cross diabatically 
the first avoided crossing aX E = 4.81kV/cm, and then 
we decrease the field (green arrows in Fig. [8]) in order to 
cross the same avoided crossing in the opposite direction 
with an intermediate velocity. Here the occupation prob- 
ability is split between the ground state (33.3%) and the 
first excited state (66.6%). This is clearly seen at time 
in Fig.[9l The electric field is further decreased, slowly at 
first and then rapidly, so that the upper branch crosses 
the complex avoided crossing aX E = and energy of 
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Figure 7: Numerical simulation of the path of Fig. [6Ua) in 
order to reach the Bell-type state {\RR) + \LL))/y/2. The 
electric field used (a) and the probabilities obtained (b) are 
shown as functions of time in semi-log plots. 



—413 meV. Then an adiabatic increase follows. We again 
cross the complex avoided crossing (time in the mid- 
dle of the avoided crossing and time after crossing it 
adiabatically) . We see that the state at is a combina- 
tion of a RL and a RR state. We continue slowly until 
the upper branch approaches the last avoided crossing at 
E — 3.02kV/cm. Here the velocity is chosen at an inter- 
mediate value, so that the occupation probability divides 
itself equally among the two states. The end result is a 
superposition of the three states mentioned at the begin- 
ning. This is clearly verified in Figs. [9] and [TOl state at 
t% = 1000 ps. In this simulation, the final square over- 
lap with the adiabatic states 1, 5, and 6 are respectively 
Pi = 0.324, Ps = 0.325 and Pg = 0.320, and then the 
final probability to find the target state is of 96.45%. 



V. FINAL REMARKS 

In this paper we extended a recently proposed con- 
trol method for quantum systems with energy spectra, 
which, as function of a control parameter, are charac- 
terized by the presence of avoided crossings. While in 
previous publications we showed how to work with di- 
abatic and adiabatic changes of the control parameter, 
that is, using the avoided crossing as a yes-no switch, 
here we explored the possibility of using intermediate ve- 
locities to cross the avoided crossings, thus obtaining lin- 
ear combinations of diabatic states. This generalization 
of our previous control strategy enables us to reach more 
general target states, not restricted to eigenstates of the 
system's Hamiltonian. 

The results presented here are for a two-electron dou- 
ble quantum-dot structure, but are not restricted to that 
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particular system. Since the issue of localization is an 
interesting one in double- well potentials, we have stud- 
ied it here as an application of our method. We showed 
how one can obtain target states with different types of 
localization, like having both electrons in one given well 
or having both electron in either well with chosen proba- 
bilities (Bell-type states) . Finally, we have shown how to 
obtain a coherent superposition of three states with each 
of the three types of localization present in this system. 
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Figure 8: (Color online) Schematic diagram of the path to 
obtain the target which is a superposition of three adiabatic 
states with different localization (Eq. |4|). For clarity, we use 
different colors when the electric field is increasing (blue and 
red) or decreasing (green). Our initial state, the ground state 
without electric field, is indicated by a filled square). The 
desired target state is a combination of the three adiabatic 
states indicated by ©. 
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Figure 9: Upper Panel: Electric field as a function of time 
used to carry out the control strategy schematically shown in 
Fig. El Lower Panel: Overlap of the evolving state with the 
first six adiabatic states. We can clearly see that the final 
state approximates well an even linear combination of the 
adiabatic states 1, 5, and 6 at the corresponding final electric 
field. 
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Figure 10: Evolving wave function at various times during the 
time evolution schematically shown in Fig. [8] and described 
quantitatively in Fig. [9l See text for details. 
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